Easy collective polarization switching in ferroelectrics 
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The actual mechanism of polarization switching in ferroelectrics remains a puzzle for many decades, 
since the usually estimated barrier for nucleation and growth is insurmountable ("paradox of the 
coercive field" ) . To analyze the mechanisms of the nucleation we consider the exactly solvable case 
of a ferroelectric film with a "dead" layer at the interface with electrodes. The classical nucleation 
is easier in this case but still impossible, since the calculated barrier is huge. We have found that 
the interaction between the nuclei is, however, long range, hence one has to study an ensemble of 
the nuclei. We show that there are the ensembles of small (embryonic) nuclei that grow without the 
barrier. We submit that the interaction between nuclei is the key point for solving the paradox. 
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The polarization switching in ferroelectrics (FEs) is 
most commonly used in applications (capacitors, memory 
elements), yet this process remains the least understood 
in spite of numerous experimental studies. As a rule, 
the ferroelectrics are switching in the field E c which is 
some order of magnitude lower than the so called "ther- 
modynamic coercive field" E cQ U-||. The latter is the 
field at which the homogeneously polarized ferroelectric 
loses its stability in external field applied in opposite di- 
rection to the polarization. The fact that the switching 
takes place well before the point of instability is reached 
means that it proceeds by inhomogeneous nucleation and 
growth of domains of a new phase. But how the domains 
nucleate? The difficulty in answering this question has 
been emphasized by Landauer in late 1950s Q]. His and 
later estimates [^]J^] showed that the energy barrier for 
creating a nucleus with reversed polarization is practi- 
cally insurmountable, U* > 10 3 fcsT in the field of about 
lOOkV/cm. This problem, or a "paradox of a coercive 
field" , has in fact been realized for domain nucleation in 
ferromagnets in 1938 ||. The experimental coercive field 
for the bulk magnetic materials was known to be many 
times smaller than that suggested by the micromagnet- 
ics theory ("Brown's paradox") The disagreement 

is especially striking for hard magnetic materials, where 
the situation is closer to the ferroelectrics considered in 
the present paper. There were numerous suggestions over 
the years that some defects can assist the nucleation and 
reduce the coercive field down to experimentally observed 
values. The situation in magnetics was summarized by 
Brown in 1965 who noted that the idea is plausible but 
"there has been no strikingly successful calculation based 
on a completely realistic model" Q]. Since then the sit- 
uation has apparently remained the same Q , whereas a 
switching in fine magnetic particles seems to be fairly 
described by micromagnetics theory. 

The goal of the present paper is to suggest a possible 



new mechanism for the polarization switching in ferro- 
electric materials. To this end we consider an exactly 
solvable model of a ferroelectric material with an ex- 
tended inhomogeneity of its dielectric response at the 
ferroelectric-electrode interface, which is called a "pas- 
sive" or "dead" layer (for references see e.g. |9"P,|lO|]). 

The main feature of the ferroelectric film with the dead 
layer is that the film exists in a polydomain state at any 
thickness of the dead layer jjj. The monodomain state 
in this system, which can be produced by cooling in the 
field, would tend to transform into the equilibrium poly- 
domain state. The transformation is favored by the fact 
that the polydomain state reduces the energy of the de- 
polarizing field in spite of increased surface energy of the 
domain walls. This is contrary to the usual notion that 
the depolarizing field hinders the nucleation Q, yet the 
classical nucleation remains impossible. 

First, we shall study the classical nucleation of the 
individual domains of the new phase. The nuclei will 
be assumed below to have a form of stripes or cylinders 
with the domain walls perpendicular to the plane of the 
film (c-domains), which is a reasonable approximation. 
We shall evaluate the barrier for their nucleation and 
show that for the individual nucleus it is practically in- 
surmountable, although the dead layer helps to reduce it. 
We then abandon the classical approach and study the 
interaction between nuclei and find that it is long range. 

It becomes clear that when an individual nucleus can- 
not grow the ensemble of nuclei may be able to. As an 
example of such an ensemble we consider a periodic ar- 
ray of nuclei and show that it indeed provides a path 
to the equilibrium state. We shall show that there is 
no energy barrier for its growth already when the nu- 
clei become larger than the thickness of the domain wall 
W. After this 'embryonic' state has been passed the free 
energy of the system decreases monotonously as the nu- 
clei grow. Obviously, the energy of embryonic nuclei is 
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much smaller than the critical energy of the Landauer's 
nucleus 0. This is simply related to the fact that the 
size of the critical Landauer needle-like nucleus is large, 
the radius of its base is r\ = WE c o/E ext ^> W, where 
E C xt R E c q) is the external field, therefore its energy 
[/£ is huge. In the present case of the dead layer E ext 
should be replaced by the (small) depolarizing field in 
the ferroelectric and the critical radius of the Landauer's 
nucleus remains ^> W. Since for the collective nucleation 
r£ oU ~ W <C r^, the incurred energy barrier for a nuclei 
to grow beyond the embryonic state £7* oll , if any, should 
be much smaller than the standard barrier [/£ for indi- 
vidual nuclei. The equilibrium density of the embryos is 
large (see below) and the nucleation of the macroscopic 
domain (and the coercive field) would be determined by 
the waiting time for optimal fluctuation and its depen- 
dence on the electric field for an ensemble which can grow 
without, or almost without, the barrier, 

We shall mainly discuss the energy aspects of the nu- 
cleation. What the present analysis demonstrates is that 
the nucleation of ferroelectric domains is facilitated by 
an interaction between the nuclei. Note that no such 
interaction is taken into account within the Kolmogorov- 
Avrami theory |0] which is widely used to treat switch- 
ing in ferroelectrics. 

We shall consider first the problem of the barrier for 
single stripe and then cylindrical domains, which appear 
to be insurmountable. We establish, however, that the 
interaction between the nuclei is long range and may give 
a clue to actual nucleation process. We then turn to 
exactly solvable case of an ensemble of stripe domains 
and show that the barrier for nucleation is actually zero 
already when the nuclei are at the embryonic stage. 

The geometry of the present problem for ferroelectric 
with dead layer is illustrated in Fig. 1. For a short- 
circuited electrodes (zero bias voltage) the free energy 
of the system is F = Fq + U es , where the electrostatic 
energy U es is 0: 



U es = \ [ dA(J(fi 
1 JFE 
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where a is the bound charge due to the spontaneous po- 
larization, ip is the electrostatic potential, while Fq in- 
cludes the surface energy of the domain walls, and the 
integration goes over the FE surface. The electrostastic 
potential if is found from solving the Poisson equation for 
assumed domain structure || . We obtain with the use of 
the Fourier transformation the total electrostatic energy 
for arbitrary one-dimensional domain structure as 
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where o~k = dx exp(—ikx)o~(x, z = 1/2) is the Fourier 
component of the surface bound charge a{x,z = 1/2), 



I the thickness of the FE film with e c ( a ) the dielectric 
constants in c (a) direction, d the thickness of the dead 
layer (Fig. 1). 

We begin with the case of a nucleus in the center of the 
plate, i.e. with the following distribution of the bound 
charge 

a(x, z = +R//2) = -(+)P S , \x\ < a/2; 

a(x, z = +RZ/2) - +(-)P s , a/2 < \x\ < R, (3) 

where 2R is the width of the plate and a is the width of 
the nucleus, and P s the spontaneous polarization. In this 
case <7fc = ^rr (— 2sin4p +sinfci?) . The dependence of 
all physical quantities on R disappears in a limit R — > oo, 
as it should. It is handy to always subtract the (constant) 
electrostatic energy of the uniformly polarized sample, 
which is characterized by the Fourier transform of the 
bound charge ou = (2P s /k) sinkR. By using Eq.([|) one 
then finds the change of the electrostatic energy due to 
creation of a stripe nucleus 
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for d < ly/e c /e a , s g = ^e a £ c - Note that this is the 
change in electrostatic energy with respect to uniformly 
polarized sample when the nucleus is present, so it does 
not apply to completely reversed sample. The electro- 
static energy favors the nucleation, since the domains 
reduce the energy of stray field . The total energy of the 
stripe nucleus per unit length is 
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where 7 = P S 2 A is the surface energy of the domain wall 
with A the temperature dependent characteristic length 
. We see that the gain in electrostatic energy eventu- 
ally overwhelms the surface energy, and there appears an 
exponentially wide barrier for the nucleus in F(a) when 
E C xt = 0. 

We shall now consider the case of two nuclei of a new 
phase in a form of stripe domains. With the help of 
Eq. (Q) one can easily calculate the change of the elec- 
trostatic energy due to formation of two nuclei having 
the same width a, with the separation r between their 
centers. We find for the energy of the interaction of two 
stripe domains per unit length for r»s 
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which corresponds to long range dipole-dipole interaction 
between two stripes |l2| . This observation indicates that 
the interactions in a system with an ensemble of nuclei 
would be very important. Note that the interaction of 



2 



the stripe nucleus with the edge of the sample (7 e dge is 
also long range, C/ c d g c(^o) = 2e~ 1 P 2 ad 2 /xq, and they are 
repelled from the edge. 

Similar treatment can be repeated for a single and a 
pair of cylindrical nuclei. The expression for the electro- 
static energy for cylindrical nucleus is similar to that for 
the stripe case (0) 
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where for the nucleus with the radius a in the center 
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of the 
radius R we obtain 07. = 2ir J Q drrJo(kr)a(r, z = 1/2) = 
(2nP s /k) [RJtQtR) - 2aJ 1 (ka)} , while for the uniformly 
polarized sampe — (2TtP s /k)RJi(kR), with J n (z) the 
Bessel function. The integral in the expression (||) for 
the electrostatic energy of the cylindrical nucleus can be 
evaluated with the result 
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The free energy of one cylindrical nucleus is 
-Fcyi(a) = 2-iral^/ + U es - 2tt P s E ext la 2 



(9) 
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It is obvious from Eqs. (|9|),(|10[) that the gain in the elec- 
trostatic energy eventually overwhelms the growth of the 
the surface energy of the domain wall with increase of 
the radius a of the nucleus. The critical radius is ex- 
ponentially large when E cyA = 0, it behaves roughly as 
a c ~ dexp (0.3a 2 K /d 2 ) , where = (0.3eA0 1/2 
Kittel period, , 
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The corresponding 



barrier height F c is very large, F c ^> P 2 d? at ~ E at , where 
d at is the characteristic "atomic" length (of the order 
of the lattice parameter) and the "atomic" energy E at 
amounts to a few eV. Thus, in this case the barrier is 
also huge, comparable to the barrier for the Landauer's 
nucleus, and its growth is prohibitively expensive. 

The interaction between two cylindrical nuclei can be 
estimated for r>oas 
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In this case it appears to also be long range of a dipole- 
dipole type 

The importance of these interactions lies in their long 
range rather than sign. Therefore, one has to accurately 
evaluate the total energy for a system of nuclei, since it 
does not reduce to a sum of asymptotic interactions (j?]) 
or (|Tl|). The long range interactions give us a clue to 
the mechanism of nucleation and growth of a new phase. 
To illustrate that it indeed may solve the problem, we 
shall consider a system of stripe-like domains. The results 



below will demonstrate that the electrostatic energy not 
only favors nucleation but eliminates an energy barrier 
for the growth of nuclei in an ensemble when their size 
is larger than the domain wall width, i.e. the switching 
proceeds collectively. 
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FIG. 1. Free energy density F/AP^ for the growth of 
ensemble of nuclei for different values L/W and the fraction 
of volume occupied by the ferroelectric f — l/L (d = L — I 
is the thickness of the dead layer). Note that there is no 
barrier for growth of very small nuclei (from 5 ~ 1) towards 
the domain pattern with equal width 8 = of the domains of 
both phases P 3 and — P a (Kittel state). The flatness of the 
free energy F in the vicinity of the Kittel state demonstrates 
remarkable softness of the ferroelectric with the dead layer. 
Inset shows schematics of the growth of the ensemble of nuclei 
(narrow domains of the new phase) . 

To study the collective nucleation we shall analyze a 
domain structure with stripe domains of opposite polar- 
ization of widths czi and 02, the period T — a\ + 02, and 
the asymmetry parameter 5 = , which measures a 

net polarization of the film ||. For zero external bias 
voltage assumed throughout the present paper <5 = in 
the equilibrium (polydomain) state, whereas in the mon- 
odomain state S = 1. At S — > 1 the system consists of 
very narrow domains with the polarization opposite to 
the net polarization (Fig. 1, inset), i.e. it is a periodic 
ensemble of nuclei. The free energy F of this system at 
zero external bias is given by 
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where D n = ^e a e c coth , 
2im/T M. Fortunately, this free energy can be found 
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analytically in the asymptotic case of narrow dead layer 
(d <C I) and correspondingly wide period of the domain 
structure, T ^> d, when e g — ^Je a e c . We have found ear- 
lier that the period of the domain structure in this case is 
exponentially large, T = 0.95dcxp(0.4a|-/<i 2 ). We note 
that in this case D n = e g (l + coth ) and summation 



in (12) can be performed to yield 
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where b = 4ird/T < 1, and Li n (z) = J2kLi z k /k n . For 
the case of an ensemble of narrow domains, 5 ~ 1 the free 
energy can be found from the known asymptotic behav- 
ior of the Li n (z) function JljJ, yielding an approximate 
expression 

F _ 2ird5 2 2AI T 2 e 3 b 2 

AP 2 ~ e c (d/l) +e g + ~ + ~ g { ~ ' n tt 2 (1 -(5) 2 - 

(14) 

Now everything depends on how the free energy F be- 
haves as a function of 5 when the nuclei grow, i.e. when 
6 reduces from unity towards zero. One can easily see 
that at x = 1 — <5 <C 1 the free energy F /AP 2 is given, 
with respect to a constant, by the function 
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where /(0) = 0. This function does not have a barrier 
as a function of x when T is kept constant. Indeed, the 
second term in (|l^) has an exponentially small maxi- 
mum (= Ts^^-Xq) at xq = Aed/T <C 1. It is, however, 
suppressed by the first term, which corresponds to the 
energy of homogeneous field created by the net polariza- 
tion and is linear in a;. As a result, there appears to be 
no barrier for the growth of nuclei. Note that we have 
actually restricted the system's path for nucleation by 
constraining the domain pattern to a fixed period. Even 
under this constraint the growth proceeds without the en- 
ergy barrier, and this would be even more so if we were 
to lift the constraint and allow the system to follow an 
optimal path to equilibrium. This behavior does not de- 
pend on the approximation we have made for evaluating 
the free energy, the exact calculation of the free energy 
( |l2| ) for all S shows that the collective growth of nuclei 
we just described proceeds without the barrier (Fig. 1). 
As mentioned above, the smallest size of nuclei where the 
present analysis applies is of the order of the domain wall 
width W and the barrier for the nucleation of such small 
embryonic nuclei is expected to be zero or much smaller 
than the usual estimates EkI for individual nucleation. 



For nucleation to proceed by the present mechanism 
the only condition is the presence of extended dielec- 
tric inhomogeneity in a sample. The external field 
does promote growth of the nuclei, but the nucleation 
in the present system occurs even without it. The 
likely requirement is that the lateral extent of this in- 
homogeneity should be much larger than the period 
of the equilibrium domain (Kittel) structure ax- One 
can estimate the rate of embryo nucleation in 1cm as 
~ (l/W 3 T p h) exp(— Uem/kT), where r p h is the charac- 
teristic (optical phonon) time, and the "atomic" esti- 
mate of the energy of the embryo is U em ~ 7VF 2 ~ 
E aty /T c /T at - 2 - 3 • 10 3 K @, with the characteristic 
"atomic" temperature T at ~ 10 4 K. Taking a conserva- 
tive estimate of the embryo lifetime as ~ 10r p ^, one finds 
the equilibrium density of the embryos ~ 10 17 cm~ 3 . At 
such high densities the embryos should 'feel' the field of 
each other, and favorable ensembles should appear within 
a reasonable time, unlike in the case of the Landauer's 
nucleus where the expectation time much exceeds the 
lifetime of the universe. 

In conclusion, we have suggested a possible way to 
solve the 11 paradox of the coercive field" by demonstrat- 
ing a collective mode of domain growth past the embry- 
onic stage with sizes about the domain wall thickness 
W, which proceeds without energy barrier. The origin 
of this cooperative phenomenon is a long-range interac- 
tion of electrostatic origin between the nuclei. The pos- 
sible screening by free charges in ferroelectric does not 
seem to be important, since the conductivity of ferro- 
electrics is usually too low to have any effect [jD. As a 
corollary, we note that the Kolmogorov-Avrami model 
(KA) jllj is inapplicable to growth of domains in ferro- 
electrics, since the essential long-range interaction be- 
tween nuclei is completely neglected in this approach 
(note that KA fully bypasses the question of how the 
domains were nucleated in the first place). The present 
results are general, and have the implication that in fer- 
roelastic materials and possibly also magnetic materials 
the nucleation would be facilitated by the long-range in- 
teraction between nuclei of a new phase. This could alle- 
viate the Brown's paradox of the coercive field in relation 
to switching in bulk ferromagnets. 
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